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The evolution of the quintessence field during a phase of chaotic inflation is studied. The inflaton <j> 
and the quintesssence field Q are described in a supergravity framework where the coupling between 
the inflaton and quintessence is induced by non-renormalisable operators suppressed by the Planck 
mass. We show that the resulting quintessence potential during inflation possesses a time-dependent 
minimum playing the role of an attractor. The presence of this attractor forces the quintessence 
field to be small during inflation. These initial conditions are such that the quintessence field is on 
tracks now. 
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I. INTRODUCTION 



A host of recent cosmoiogicai observations, the 
anisotropy of the Cosmic Microwave Background 
(CMB) yl 0, the large scale structures of the Uni- 
verse Q and type la supernovae Q, indicate that the 
Universe has experienced two stages of cosmic acceler- 
ation. The first one is the inflationary era which oc- 
curred in the early Universe It is responsible for 
the almost flatness of the Universe and primordial den- 
sity fluctuations @ (see also Refs. 0). The second one, 
which started in the recent past, leads to the present 
acceleration of the expansion of the Universe. Various 
explanations for this last phenomenon have been pro- 
posed in the literature: a pure cosmoiogicai constant 

quintessence n m m Gun q m, k-essence lb , 

modified gravity theories \u\. the Chaplygin gas HS, 
bulk viscosity or quantum cosmoiogicai effects |20j . 
In this paper, we focus on the quintessence hypothesis. In 
this case, the two phenomena described above are mod- 
eled as resulting from the presence of two scalar fields 
whose energy densities drive the acceleration of the ex- 
pansion. The quintessence hypothesis has been further 
investigated in Refs. |2l|. In particular, finding a natu- 
ral candidate for the quintessence field in the realm of 
high ener gy phys ics has been a major goal for lot of au- 
thors 0,13 as wen as studying some aspects of 
its interaction with the "rest of the world" [23]. Since, 
contrary to a cosmoiogicai constant, the quintessence 
field can develop some inhomogeneities, the theory of 
cosmoio gica i perturbations has also been studied in de- 
tails (l4[l24j and has been used in order to constraint var- 
ious models observationally [2j|. Of course, the prospect 
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of utilizing the fact that the quintessence equation of 
state is no longer time (or redshift)-independent as a tool 
for discriminating amongst the various possibilities has 
been widely discussed in the recent literature 26] . 

As a low energy description of string theory |23|, su- 
pergravity captures prominent features of physics beyond 
the standard models of particle physics and cosmology. 
Supergravity is the best framework within which both 
quintessence and inflation can be described. Indeed in- 
flation (in its most common models like chaotic infla- 
tion) involves high energies as the inflaton rolls down 
its potential with values exceeding the Planck mass. 
Similarly, in quintessence models with a rolling scalar, 
the quintessence field reaches values of the order of the 
Planck mass now. Hence the necessity for a treatment 
where non-renormalizable interaction terms suppressed 
by the Planck mass are under control. In supergravity, 
such non-renormalizable corrections to supersymmetric 
models are taken into account and play an important role. 
This justifies the use of supergravity models both in in- 
flation and quintessence model building. In the following 
we will concentrate on both quintessence and inflation as 
described in supergravity. 

One of the commonly used models of quintessence, first 
devised by Ratra and Peebles requires an inverse 
power law behavior V(Q) = M 4+a Q~ a with an attrac- 
tor mechanism at large time. It was soon realized that 
this type of potential can be generated in supersymmetric 
theories when a strongly interacting sector is present |ll| . 
In particular, the value of the quintessence field becomes 
of the order of the Planck scale which prompts the ne- 
cessity of a supergravity treatment. A simple embedding 
of the previous model in supergravity fails as the poten- 
tial is highly modified by supergravity corrections and 
can become negative 0, |p| . Hence a more phenomeno- 
logical approach may be required where one postulates 
the form of the Kahler potential and the superpotential 
which leads to quintessence in supergravity. This was 
done in Refs. 0,0] and subsequent work. 

Similarly, as a high energy phenomenon occurring in 
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the early universe, inflation must be described within su- 
pergravity. Recently, there has been a upsurge of in- 
flation models in supergravity motivated by string the- 
ory [2ll29ll30|. It seems natural to study the influence 
of quintessence on inflation and vice versa. 

Quintessence must be almost decoupled from ordinary 
matter, otherwise the quintessence field would lead to 
observable fifth force signals |23j. On the contrary the 
inflaton field must couple quite strongly to ordinary mat- 
ter in order to have a reheating period at the end of 
inflation where the oscillations of the inflaton result in 
a radiation bath. Hence the coupling of quintessence 
and inflaton cannot be large. A natural way of real- 
izing this criterion is to consider that the inflaton and 
the quintessence field arc decoupled in the Kahler poten- 
tial and the superpotential of supergravity. This implies 
that the only possible interactions between both fields 
spring from non-renormalizable interactions suppressed 
by the Planck mass. Here we provide such a super- 
gravity description of the coupling between inflation and 
quintessence. 

The paper is arranged as follows. In a first part 
(Sec. II), we analyze the supergravity coupling between 
a particular quintessence model, the so-called SUGRA 
model [T2I fl3l Il5| , and a generic inflationary model. We 
then (Sec. Ill) apply this analysis to the specific example 
of chaotic inflation where we show that the quintessence 
field develops a potential with a rolling minimum during 
inflation. The rolling minimum is an attractor such that 
the values of the quintessence field remain small through- 
out the inflationary era. In particular, these values are 
much smaller than the values of a free quintessence field 
during inflation. The smallness of the quintessence field 
during inflation implies that it is on tracks now, i.e. it 
reaches its long time attractor. Finally, in Sec. IV, we 
discuss the limitations of our approach, try to indicate 
what possible improvements could be and present our 
conclusions. 



It is worth mentioning that one can derive the SUGRA 
model from more general Kahler potentials but we will 
not need them in this article, see Ref. [l5|. At this stage, 
we assume that 



(X} = £, (Y}=0. 



(3) 



As a specific example, £ can be realized as a 
Fayet-Iloupoulous term arising from the Green-Schwarz 
anomaly cancellation mechanism. When /i = 0, the Q 
direction is flat. It is lifted by the superpotential leading 
to the quintessence potential. In supergravity, negative 
contributions to the scalar potential arise from the vac- 
uum expectation value (vev) of the superpotential. No- 
tice that we have here 



quint/ 



o. 



(4) 



implying that no negative contribution appears in the 
scalar potential. 

We are now in a position where the scalar potential 
can be computed. In supergravity, it is given by 



V 



1 



e G (G A G A - 3) 



(5) 



where the matrix G A g which is used to raise and lower 



the index A is defined by 

d 2 



G 



AB 



d<p A d(tp B ) f 



KK quint + In 



(- 3 | 



TI'« 



quint | 



(6) 



where tp A = {X,Y,Q} are the fields in the quintessence 
sector. Straightforward calculations leads to a matrix 
which is block-diagonal, namely 



^ab — K 



1 

(kQQ^) p 
1 



(7) 



Then the complete SUGRA potential becomes 



II. QUINTESSENCE AND INFLATION IN 
SUPERGRAVITY 

A. Quintessence in Supergravity 

Let us now briefly review a simple model of 
quintessence in supergravity often dubbed the SUGRA 
model in the literature. We assume that the Kahler po- 
tential and the superpotential are given by 0, ^| 

Squint (X, Y, Q) = AAt + QQt + k p YY \ (QQ\y (i) 



quint 



(X,Y,Q) = f,X 2 Y. 



(2) 



with k = 8ir/m 2 r Here X and Y are two charged fields 
under an (anomalous) U(l) symmetry with charges 1 and 
—2, while Q is the neutral quintessence field. Notice the 
direct coupling between Q and Y. The constant /i is a 
dimensionless coupling constant and p is a free coefficient. 



V^juint(Q) 



M 4+2p 



(8) 



where the mass scale M characterizing the potential can 
be expressed as M 4+2p = h 2 £ i a k~ p . It is easy to find that 
(-Kquint)gQt = 1 which means that the real part field Q 
is in fact not correctly normalized. Therefore, one has 
to redefine the field Q according to Q — > Q/y/2 and this 
gives 



U quint (Q)=e^ 2 / 2+ ^ 



M i+2p 

~Q 2p ~ 



(9) 



where we have slightly redefined the mass scale M such 
that M 4+2p -> M 4+2p x 2P. The main feature of the 
above potential is that supergravity corrections have 
been exponentiated and appear in the prefactor. Phe- 
nomenologically, this potential has the nice feature that 
the equation of state u> = Pq/ Pq can be closer to —1 than 
with the Ratra-Peebles potential. 
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B. Inflation in Supergravity 



case, the potential takes the form 



Let us now give a brief description of the inflation mod- 
els we will concentrate on. We will consider a class of 
models described by the following Kahler potential 



Ks 



inf 



= In 

K 



K 1 ' 2 (p + p*) -KK(<f>- tf) 



(10) 



where K, and Q are arbitrary functions. The field <p 
is the inflaton while p represents, for instance, a mod- 
uli of a string compactification. The superpotential 
Wi n f = Winf {p, 4>) is not specified at this stage. This form 
is justified by the fact that one can obtain flat enough po- 
tentials in supergravity by requiring that a shift symme- 
try <j> — ► <j) + c, where c is a real constant, is a symmetry 
of the Kahler potential, later broken mildly. The Kahler 
potential given by Eq. (|10|l obviously possesses this sym- 
metry. Indeed, a striking feature of F-term inflation in 
supergravity is the natural presence of 0{H ln {) correc- 
tions to the inflaton mass which would spoil the flatness 
of the potential. These problems can be avoided by con- 
sidering the above class of models. 

To go further, one must specify the functions /C, Q and 
the superpotential. We choose an example of chaotic 
inflation as can be found in Ref. |3ll | where a similar case 
is treated. Explicitly, one assumes 



and for the superpotential 



+ \{<j>-^)\ (11) 



a 

—mo 

2 V 



(12) 



The factor a in the superpotential is free and can be 
chosen for future convenience. Notice that the shift sym- 
metry is preserved by our choice of the functions K. and 
Q while, on the contrary, the superpotential breaks this 
symmetry explicitly. Then, straightforward calculations 
lead to 



1 



A 2 (3 



2 2 ±2 

-a m ((> 



(13) 



where A = K 1 / 2 (p + pf). It is easy to see that the mod- 
uli can be stabilized if A = 2. Furthermore, one can 
check that the normalization of the inflaton is given by 
(Kmf)^ = 3/A — 1 = 1/2 and, hence, is correct. In this 



V ini {4>) = T m 2 4> 2 , 



(14) 



which is nothing but the usual chaotic inflation potential 
if one chooses a = 

Of course, it is possible to discuss more complicated 
and/or general inflationary models. The one considered 
here has the advantage to lead to the prototypical single 
field inflationary model, namely chaotic inflation. Since 
our main goal is not to study inflation itself but the 
coupling of the inflaton with the quintessence field, this 
model is sufficient. However, it is clear that the next step 
would be to study how the form of the coupling term that 
we are going to derive below depends on the assumed in- 
flationary model. 



C. Coupling the Inflaton to the Quintessence Field 

We now turn to our main goal, namely the calcula- 
tion of the coupling between the inflaton field and the 
quintessence field. Our basic assumption is that the 
quintessence and inflation sectors are decoupled, i.e. that 
the total Kahler potential and superpotential can be writ- 
ten as 

K = K quint {X, Y, Q) + K in{ (p, oh) , (15) 
W = VF quint (X, Y, Q) + W in f (p, <f>) , (16) 

where the quintessential Kahler potential and superpo- 
tential have been given before but where, at least at this 
stage, the inflationary part is still arbitrary. However, 
later , we will restrict our considerations to the (chaotic) 
inflation model studied in the preceding section. From 
the above equations, one deduces that the matrix G A §, 
where now ip A = {X,Y,Q, p,<j)}, is diagonal by blocks. 
Explicitly, one has 





Gppi 










0" 

















Gab — 








K 
















K 





















K_ 



(17) 



Then, the scalar potential takes the form (recall that the 
D-terms contribution vanishes) 



V 



^ Q2/2 v inl (p,, 



V quint (Q) + K 2 [Z 2 
I 



9Ly Wint \ 2 e< K ^ 2 / 2 ) 



(18) 



where, 



Vi a { (p, (f>) 

Kjuint (<3) 



3 G lnf 



7 (— [G£,(Gtaf) A -3] , (19) 



Let us notice that, for convenience, we have slightly 
changed the notation for V qu i nt (Q). Now, we no longer in- 



_ n «Q 2 /2 



M 4+2 p 



(20) 
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elude the factor exp(/-t£ 2 ) in its definition, see Eq. JSJ. As 
explained before, we have also redefined the quintessence 
field according to Q — > Q/V% in order to work with cor- 
rectly normalized fields. The above expression represents 
the general form of the coupling between the SUGRA 
model of quintessence and inflation in supergravity. 

We now specify the inflaton model and consider the 
model described in the previous subsection with a = y2- 



It is convenient to work in terms of dimensionless quanti- 
ties. In particular, we define the dimensionless potential 
/(</), Q) by V((f>, Q) = TOp i /(</',Q) and this quantity can 
be written as 

/ {4>, Q) = fixd + /quint + /inter , (21) 

where 



/inf = 



1 / m 



2 m„ 



4:TT Z 



C W PI V 7 Pi , / quint 



1 ( M \ 

■i / 



i+2p 



8 V ITl- 



_Q_ 



-2p 



1 (JL 

2 \ m„, 



(22) 



At this point, some remarks are in order. A priori, there is no clear separation between the inflaton and the quintessence 
fields in the term f- m { because of the presence of the exponential term. However, in the regime we will be studying 
(during inflation), Q <C m Pl and, therefore, the exponential term will be very close to one. In this case, one recovers the 
simple chaotic model V m f = m 2 </> 2 /2. The term / qu int is nothing but the SUGRA potential studied in Ref. [T2I IT2I ITsT ] 
but, during inflation, it will reduce to the Ratra- Peebles case. Let us notice that we have an extra factor 1/8 originating 
from the term exp (kK^). This comes from the fact that nK ln { = — 31n(2) since the moduli is stabilized at A = 2. 
Finally, in the regime Q <C m pl , the interaction term reads Vintcr cx rn^^Q 2 /m^ r This is due to the fact that we 
have Q 3> £ as will be discussed below. The coupling constant between the inflaton and the quintessence fields reads 
to V to p ; - The Planck mass appears in this expression because the coupling between the two fields has been entirely 
fixed by the supergravity. Notice also that the quintessence field picks up an inflaton dependent mass term during 
inflation. The competition between this mass term and the Ratra-Peebles potential will be studied in the next section. 
Finally, using the fact that 



Of 



d((f>/m Pl ) 

Of 
d(Q/m Pl ) 




(23) 



it is easy to show that this potential possesses an absolute 
minimum given by 



III. COSMOLOGICAL EVOLUTION 



A. Fixing the free parameters 



= 0, 



_Q_ 

m D1 



(25) 



Let us now discuss the values of the free parameters 
that appear in Eqs. (|21|l and l|22|l . If we assume that the 
quintessential part is responsible for the acceleration now 
then one should have 



The potential is represented in Fig. ^ The value <fi = 
is the minimum of the inflaton potential without interac- 
tion while Q = \Jp/{4,-k) is the minimum of the SUGRA 
potential. At the absolute minimum, the value of the 
potential is non-vanishing and given by V — Vq U i nt = 
(m* i /8)(M/m pl ) i+2 P(4ir/p)P exp« 2 ) exp(p). 



M 4+2p 



(26) 



where Qo and Hq denote the values of the quintessence 
field and of the Hubble parameter (respectively) now, at 
vanishing redshift. In order to have a successful model of 
quintessence, the field should be on track today which in 



5 




FIG. 1: Upper panel: Potential V(cf>,Q) for the following choice of parameters: p — 3, m = 10 5 m pl , £ = 10 30 m pi , and 
(M/m P! ) 4+2p = 10~ 122 . The absolute minimum located at Q ~ 0.4886 x m pi and = cannot be viewed with the scales used. 
Bottom panel: zoom in the region of the potential where the minimum is located. It is clear that the tiny values of the inflaton 
field are, in this panel, not interesting from a physical point of view (there is no inflation for such small values). 
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turn implies that Qq = 0(m pl ). Strictly speaking, this 
conclusion is valid for the Ratra-Peebles potential only, 
but it has been shown in Ref. 0, ^1 that this is also 
valid for the SUGRA potential despite the presence of 
the exponential correction (this is simply because, except 
at small redshifts, we have Q <C m pl and the exponen- 
tial SUGRA correction does not play an important role). 
This gives 



M_\ 4+2p 
m Pi) 



Hi 



10 



-122 



(27) 



Using the fact that M A+2p ~ /i 2 £ 4 ft _p and assuming no 
fine-tuning of the coupling constant, i.e. /i = 0(1), one 
deduces that 




(28) 



In a sense this is the usual fine-tuning of the cosmological 
constant, it reappears here in the guise of the tuning of 
the vev of a field leading to the quintessence potential. 
However, if one works with an effective model valid u 
to a cut-off scale m c , then it has been shown in Ref. 
that the previous problem can be solved provided the 
scale is chosen such that m c <C m pl . Then the value 
of the Fayet-Iliopoulos term can even be above the weak 
scale. However, again, our purpose here is not to study 
the details of the dark energy model and, therefore, in 
the following, we will ignore these subtleties and work 
with the value of £ derived before. Let us notice that 
even with a Fayet-Iliopoulos term above the weak scale, 
in general, we still have Q > £ and then the form of the 
coupling is not modified when one works with a cut-off 
scale much below the Planck scale, see also the discussion 
after Eq. 

We now discuss the constraint on the parameter char- 
acterizing the inflaton sector, i.e. the mass to of the field. 
In order to simplify the discussion, we will assume that 
the initial conditions are such that quintessence field is 
always subdominant. In this situation the quantum fluc- 
tuations of the inflaton field are at the origin of the CMB 
anisotropy observed today. As a consequence, and as is 
well-known, the COBE and WMAP normalizations fix 
the coupling constant of the inflaton potential, namely 
the mass m in the present context. More precisely, for 
small £, the multipole moments are given by 



C. 



2H? 



1 



25em 2 £{£ + 1) 



(29) 



and what has been actually measured by the COBE 
and WMAP satellites is Q 2 ms _ PS /T 2 = 5C 2 /(4tt) ~ 

(18 x 10" 6 /2.7) 2 ~ 36 x KT 12 . The quantity H in{ is 
the Hubble parameter during inflation and is related to 
the potential by the slow-roll equation H? n{ ~ /«V[nf/3 
evaluated at Hubble radius crossing. Putting everything 
together, we find that the inflaton mass is given by 



45tt N.. 



Q 



rms-PS 
~J^2 



(30) 



that is to say 



1.3 x 10" 



X TO, 



(31) 



All the parameters of the potential are now specified. 

Let us now discuss in more detail what are the condi- 
tions under which the inflaton field is always dominant. 
The quintessence energy density must be smaller than 
the inflaton energy density. This gives a lower bound on 
the possible values of the field Q which can be expressed 
as 



Q 



low 



m <f) 

Hq TO pi 



-1/p 



1Q-66/P ( J_ 



-i/p 



, (32) 



where we have used the value of m ~ 10~ 6 x to pi ob- 
tained before. The fact that we obtain a lower bound is 
consistent with the fact that the potential is an inverse 
power-law of the quintessence field: the smaller the field 
is, the larger the corresponding energy density is. Sec- 
ondly, there exists also an upper bound coming from the 
fact that the interaction energy density must be smaller 
than the inflaton energy density. Concretely, this gives 



Qur 




(33) 



where we have used the fact that the maximal value of 
the inflaton field is <f> ~ 10 6 x to P! , see below. 

Under the condition that Qi ow < Q < Q up , the be- 
havior of the background is determined by the energy 
density of the inflaton and it is well-known that, in this 
case, the slow-roll approximation is valid. The slow-roll 
approximation is controlled by two parameters (in fact, 
at leading order, there are three relevant slow-roll param- 
eters but we will not need the third one) defined by 



H 



2He 



(34) 



In the present context, where the inflaton potential is 
proportional to m 2 (j) 2 , the slow-roll parameters are given 
by 



1 



2JV, 



1 



6 = 0. 



(35) 



where N* ~ 60 is the number of e-folds between the 
time at which scales of astrophysical interest today left 
the Hubble radius during inflation and the end of infla- 
tion. In the situation where these parameters are small, 
namely £ < 1 and 5^1, the equation of motion of the 
inflaton field can be easily integrated. For this purpose, it 
is convenient to express everything in terms of the num- 
ber of e-folds (not to be confused with TV*) defined by 



N = In 



Oil 



(36) 
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such that, at the beginning of inflation, one has N = 0. 
Then, in the slow-roll approximation, one obtains that 
the evolution of the field is given by 



(37) 



where (/>;„; is the initial value of the field. This value is 
related to the total number of e-folds given by 




(38) 



If N m in is the minimum number of e-folds required in 
order to solve the problems of the hot big-bang model 
(-Wmin — 60) then one has 



0ini > m 




3.1 x m. 



(39) 



There exists also an upper bound for the value of the 
inflaton field which corresponds to the situation where 
the potential energy density m 2 2 /2 is Planckian. Using 
that m ~ 10~ 6 x m pi this immediately gives that max ~ 
10 6 x m pi . This is the value of </> max that we considered 
before. 



B. Analytical Study of the Klein-Gordon equation 

We now turn to the resolution of the quintessence equa- 
tion of motion, i.e. the Klein-Gordon equation. It can be 
written as 



Q + ZHQ+-^V{^Q) = 0. 



(40) 



with H = a I a is the Hubble parameter which only de- 
pends on the inflaton energy density. A dot denotes a 
derivative with respect to cosmic time. We now work 
with the new time variable introduced before, namely 
the number of e-folds. One gets 



_Q 

dN 2 [m r 

"if J d(Q/m rl ) 



1 AH 



d 

In 



of 



(41) 



The main feature of the above equation is that the 
potential is now explicitly time-dependent because of 
the interaction of Q with the inflaton, i.e. we have 
/ ((/>, Q) = f (N, Q). This renders this equation very dif- 
ficult to solve exactly. Therefore, in order to get some 
analytical approximate solution, it is necessary to make 
some assumptions. With a very good accuracy, the fact 
that Qiow < Q < Q U p implies that the potential con- 
tains two dominant terms and can be approximated as 
~ Q~ 2p + Q 2 ■ As a consequence, its derivative can ex- 
pressed as 



Of 



p ( M \ 4+2p / Q 



d(Q/m Pl ) 



4 V m-i: 



-2p-~L 



+4tt 2 



Q 



m„, / m p , \m pl 



(42) 



This means that the quintessence field evolves in a time- 
dependent potential which possesses a minimum. Of 
course, this minimum is explicitly time-dependent and 
can be expressed as 



Qmm(A0 = m p( x 



16tH 



As shown below, this equation turns out to be one of 
the main result of the present article. Indeed, we will 
demonstrate that, after a period of rapid oscillations, 
the quintessence field always tends toward the above so- 
lution. Therefore, in the case where the interaction be- 
tween the inflaton and the quintessence field is important, 
Qmin can be viewed as a kind of attractor solution since, 
regardless of the initial conditions, the final value of the 
field is always given by Q min . 

Several remarks are in order at this stage. Firstly, 
let us evaluate the typical time of evolution of the min- 
imum. It is given by AN ~ Qmin/(dQ m i„/dA r ) = 
[dlnQmin/dA/] -1 . From Eq. O, one has Q min cx 



rn 
m 



¥■1 



<P(N) 



l/[2(p+l)] 



(43) 



I 

#-2/(p+i). Therefore, this implies that 

p+1 



AJV„ 



2c 



> 1 



(44) 



where e = — (dH/dN)/H is the first slow- roll parameter. 

Secondly, it is interesting to calculate the effective mass 
of the quintessence field at the minimum of its time- 
dependent potential. Using Eq. 1)43)1 . one obtains 



'off 



dp 



m 2 d{Q/m rl ) 2 



( m \ 2 


U(N)1 


\ m Pi) 





(45) 
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Therefore, one has 



functions |32| 1 and the quantity x is defined by 



2/3 



<KN) 



(54) 



'cir 



H 2 



> 1, 



(46) 



and we conclude that, at its minimum, the quintessence 
field is not a light field. 



Initially, and during a few e- foldings, one has i» 1. In 
this case, one can use the asymptotic expression of the 
Airy function and one obtains 



SQ 



e -3Ar/2 7r -l/2 a .-l/4 



Ax sin ( 7^x 3/2 + j 



+A 2 cos ( -x 



2 3/2 , * 



(55) 



From this expression, one sees that one has damped oscil- 
lations. The period of the oscillations can be very easily 
estimated. One has 



Thirdly, we are now in a position where one can 
study how small fluctuations behave around the time- 
dependent minimum. The fluctuations are given by 



SQ = Q 
equation 



Qmi n and their evolution is governed by the 



diV 2 



SQ 

TO , 



dN 



SQ 



meg 
H 



SQ 



0,(47) 



where, in the damping term, we have neglected the 
derivative of the Hubble parameter which is nothing but 
the slow-roll parameter e. Using the expression of the ef- 
fective mass established before and the expression of the 
inflaton in the slow-roll approximation, one finds that the 
solution can be expressed as 



SQ 



e -3A72 [^ 1 Ai(-x)+A 2 Bi(-x)], 



(48) 



where A± and A 2 are two constants determined by the 
initial conditions. The functions Ai and Bi are the Airy 



AJV ra 



2tt 



v/3(p+T) T 



AT-i/2 



(56) 



1 Another method to solve the Klein-Gordon equation, under the 
assumption that the slow-roll hypothesis is valid for the inflaton, 
is the following. Instead of directly neglecting e = — (dH/dN)/H 
in the damping term of Eq. 1411 as we did before, one works with 
the dimensionless field q(N) defined by 



g(N)g(N) 



H 



-1/2 



-3N/2 



q(N). 



(49) 



Then, from Eq. 1411 . it is easy to show that the field q(N) obeys 



d 2 q 
dN 2 ' 

1 



1 d 2 H 1 f dH \ 2 15 dH 9 
2H dN 2 + AH 2 \dNJ ~ THdN ~ 4 



9{N) 



H 



Of 



d(Q/m pl } 



0. 



Q/m =g{N)q(N) 



(50) 



In the second term between squared brackets, the various deriva- 
tives of the Hubble parameter can be expressed in terms of the 
slow-roll parameters. Explicitly, this term reads e(3e + 2<5)/2 + 
e 2 /4 + 15e/4 - 9/4 ~ -9/4. Therefore, the Klein-Gordon equa- 
tion can be simplified further and we obtain 



d 2 q 
dN 2 



9 1 

-o H 

4 q g(N) 



H 



df 



d(Q/m pl ) 



g(N)q(N) 



:0. (51) 



Using the fact that the potential is given by m 2 g<f> 2 /2, hence we 
are now studying SQ = g(N)q(N), the above equation takes the 
form 



d 2 q 
dx 2 



+ xq = 0, 



(52) 



where x is defined in Eq. 1541 . As before, this equation can be 
solved in terms of Airy functions and this gives 



H 



-1/2 



-3N/2 [BiM { __ x) + B2Bi { _ x)] _ (53) 



This solution should be compared with Eq. 1481 . We see that the 
equation are similar up to the factor (H/m rl )~ 1 / 2 . As our ap- 
proximation is valid during a few e-folds only, the Hubble param- 
eter can be considered as a constant and then the two solutions 
are identical. 
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where we recall that N T is the total number of e-folds 
during inflation. For a typical model with p — 3 and 
N T ~ 60, one gets AA osc i ~ 0.24. The previous equation 
also means that if the inflaton field starts at large values, 
then the period of the oscillations will be extremely rapid. 
For instance, if inflation starts at Planckian density, then 
(j>i n i ~ 10 6 m p; which implies that N T ~ 10 12 . As a con- 
sequence, one can get values as small as AiV osc i ~ I0~ 6 . 

Therefore, from the above considerations, one reaches 
the conclusion that 



AN U 



AN n! 



= 0(1)V_JL » i 



(57) 



This means that the oscillatory phase is very quick in 
comparison with the typical time scale of evolution of 
the minimum. To put it differently, the minimum can 
be considered as motionless or as "adiabatic" as the field 
rapidly oscillates and quickly joins its minimum. There- 
fore, Qmin(N) can be viewed as an attractor since it does 
not depend on the initial conditions for the quintessence 
field. 

We have established the above result under the as- 
sumption that the initial deviation from the attractor 
Q m i n is not too large (or, in other words, that SQ is not 
too large). What happens if this is not the case, i.e. if 
IQini — Qmin | > 1 or <5<3i n i 3> 1? Is the attractor still 
joined rapidly enough (i.e. before the end of inflation)? 
A priori, to answer this question requires a full integra- 
tion of the equation of motion (or a numerical integra- 
tion, see the next subsection) which is not possible. How- 
ever, we can gain some partial insights using the follow- 
ing considerations. If one has -C Q m in then the term 
proportional to Q~ 2p dominates in the potential and the 
Klein-Gordon equation remains non linear hence difficult 



to integrate. But if we now assume that we start from a 
situation where Qi n i ^> Q min , then the term proportional 
to Q 2 dominates in the potential. As a consequence, the 
derivative of the potential can be written as, see Eq. H42|) 



Of 



d(Q/m m ) 



4:TT Z 



£L) '(-*-)'*-. (58, 



Therefore, the Klein-Gordon equation is now linear and 
can be integrated. In fact, one obtains the same potential 
as before for 5Q, up to an unimportant factor 2(p + I), 
except that now Q — Q mm > needs not to be small. 
As a consequence the solution will read the same, that 
is to say, roughly speaking, Q ~ Qini exp(— 3N/2), this 
solution being valid provided Q ^> Qmin- 

Equipped with this solution, one can now estimate how 
many e-folds are necessary for the field to roll down the 
potential from a given initial condition and to reach the 
region of the minimum. When the field enters this re- 
gion, the previous solution is no longer valid because the 
term Q~ 2p starts playing a role but, on the other hand, 
since the field is now close the Q m m the calculation of SQ 
applies. In order to get an upper bound on the number 
c-folds, let us assume that Q is initially as far as pos- 



sible from the minimum, i.e. Q Y 



Q u 



pi 



Eq. (El . 



Then the number of e-folds N is solution of the algebraic 
equation 



-3AN/2 



(59) 



The solution can be exp ressed in terms of the 
Lambert function Wq |33j defined by the relation 
W(z) exp[W(z)] = z. Explicitly, one obtains 



J 



AN = 



2 fa / 6- ■ \ ( 3+ p)/( 1+ f) 

_(p + l)tf T + JWb | (p+ i)iV T (jgij 



p 



167T 2 TO 2 



-l/[2(p+l)] 



3 3(p+l)iV T /2 



(60) 



r 



Since the argument of the Lambert function Wq is very 
large, one can use the approximation Wo(z) ~ ln(z). In 
this case, one gets 



AN ~ - In 
3 



2 ( 3 +P) ln ^ 



3(p+i; 



in 



3(p + l) 



3(1 +p) 
m 
m 



pi 



In — . (61) 



For the fiducial model with p = 3, (j)- m i — 3.1m pi and 
N T = 60 one obtains AN ~ 26 < 60. Therefore, even 
in the extreme case where the quintessence field starts at 
Q up , inflation lasts enough e-folds so that Q has time to 
reach the attractor. 



In conclusion, in this subsection, we have shown that, 
during inflation, the evolution of the quintessence field is 
characterized by two very different time scales. One scale 
describes the evolution of the adiabatic time-dependent 
minimum while the second one represents the period of 
the rapid oscillations around this minimum. We have 
demonstrated that Q m in is in fact an attractor and that, 
regardless of the initial conditions, the quintessence field 
always has enough e-folds during inflation to join this 
attractor. Although the above conclusion has been es- 
tablished in the quadratic part of the potential, it is in 
fact true even in the regime where the potential is pro- 
portional to Q~ 2p as confirmed by a numerical study of 
the Klein-Gordon equation. 
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C. Numerical Study of the Klein-Gordon Equation 

We have just seen that the equation of motion cannot 
be analytically integrated with the complete potential. 
As a consequence, the part where the potential is pro- 
portional to Q~ 2p has not been explored for values of the 
initial conditions far from the minimum. In this subsec- 
tion, we perform the integration numerically. The diffi- 
culty is that we have to deal with very small quantities. 
It is therefore necessary to absorb these small quantities 
into a redefinition of the quintessence field which greatly 
facilitates the numerical integration. For this purpose, 
we write 



Q 



(62) 



where A is a constant. It is easy to show that, if A is 
chosen to be 



A 



in 
m 



pi 



-i/(p+i) 



i/(p+i) 



(63) 



then we can remove the dangerous coefficients from the 
equation of motion which now reads 



d 2 Q 
diV 2 



3 IdiA dQ A 

HdN J dN 4tt V m F 



4 



= 0. 



(64) 



In particular, it is interesting to evaluate how the time- 
dependent minimum looks like after the rescaling. From 
Eq. (I4l5|) . one gets 



Q^(N) = (^) 



l/(2p+2) 



-2/(p+l) 



(65) 



The results of the numerical integration are presented 
in Figs. El and |21 for two different initial conditions. We 
always assume that the inflaton field starts at </>;„; = 
3.1m pi which, as already mentioned, means that the 
total number of e- folds during inflation is N T = 60. 
The quintessence potential has been chosen such that 
p = 3. This also completely specifies <2 m i n and in partic- 
ular we have Q m in{N — 0) ~ 0.34 as can be checked 
directly on the figures. The evolution of Q m i n (N) is 
represented by the red dotted curve in Figs. [3 and 
Moreover, with the values of Hq and m discussed be- 
fore, the rescaling constant A is equal to A ~ 1.6 x 10~ 14 
and, therefore, the initial value of the attractor is in fact 

Qmin{N = 0) ~ 5.6 X 10- 15 - 



In Fig. [3 one has Q m i — 0.05 or Qi, 



8.3xl0- 16 m 



Initially, the field is therefore in the region where the po- 
tential is proportional to Q~ 2p , i.e. the region which was 
explored analytically before. We see that the evolution 
is very similar to what was discussed before. We have a 
period of rapid oscillations and then, after a few e- folds, 



the attractor is joined. We notice that the period of 
these oscillations is in full agreement with the estimate 
of Eq. I|56|l . One can check that the amplitude of the 
oscillations decreases as exp(— 3N/2) as demonstrated in 
the previous subsection. We conclude that all the prop- 
erties established before are confirmed by the numerical 
study, even in the part of the potential where it is pro- 
portional to Q~ 2p . However, one should also notice that, 
if the field is initially very displaced from its minimum 
such that Qini <C Q m in, then the simple Fortran code 
used to integrate the equation of motion can quickly run 
into numerical problems. This is probably due to the fact 
that Q~ 2p is a very steep potential. Despite this remark, 
one sees no reason why, in this regime, the evolution of Q 
should be different from what has been described before. 



In Fig. El one has Q ini = 20 or Q ini ~ 3.3 x 10" 
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This time, one starts from the other part of the potential, 
where Q ln > Q m in an d V(Q) oc Q 2 . The remarks made 
before also apply to this case which appears to be in full 
agreement with the analytical estimates of the previous 
subsection. 



IV. DISCUSSION AND CONCLUSIONS 



In order to study the influence of the interaction term 
and to compare its effect with the standard case, it is 
interesting to give the evolution of the quintessence field 
when this one does not interact with the inflaton (and 
when Q remains a test field). In particular, we are inter- 
ested in calculating, for a given initial condition at the 
beginning of inflation, the value of Q at the end of infla- 
tion (or at the beginning of the radiation dominated era) 
in both cases (i.e. with and without interaction). The 
case without interaction can be easily treated because the 
Klein-Gordon equation can be integrated in the slow-roll 
approximation, see Ref. |34|. In fact, this equation can 
be re-written as 



Q 



H{<f>)Q 3iP(0) 



(66) 



where a dot denotes a derivative with respect to cos- 
mic time. In the above equation V^p (Q) now means the 
Ratra-Peebles potential, namely V RP (Q) oc Q~ 2p since 
this is the potential for the quintessence field in absence 
of any interaction with the inflaton field. Due to the 
smallness of the parameter e, the slow roll approximation 
can be applied to the equations describing the motion of 
the quintessence field without interaction if the following 
condition is satisfied 



C(Q) 

3H 2 {(/)) 



« 1. 



(67) 
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FIG. 2: Upper panel: Evolution of the quintessence field during inflation (solid black line). The model of inflation is chaotic 
inflation with a massive potential and the initial value of the inflaton is chosen to be <pi n i = 3.1 x m pi corresponding to a total 
number of e-folds N T — 60. The potential of the quintessence field is of the Ratra-Peebles type with p = 3. The initial value 
of the quintessence field is taken to be Qi n i = 0.05 or Qi n i ~ 8.3 x 10 _16 m pl . This initial value is such that Qj n j < Q m i n where 
Qmin is the time-dependent minimum of the effective potential. The evolution of Q ln i n (N) is given by the dotted red curve. 
Bottom panel: a zoom of the upper figure at the beginning of inflation. 
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FIG. 3: Upper panel: Evolution of the quintessence field during inflation (solid black line) with the initial condition Q ln i — 20 or 
Qini — 3.3 x 10 _13 m P! . This initial condition corresponds to a situation where Qmi > Qmm. The parameters characterizing the 
model are identical to those used in Fig. [5] The dotted red curve represents the time-dependent minimum Qmin(N). Bottom 
panel: a zoom of the upper figure at the beginning of inflation. 
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If one applies this condition to the Ratra-Peebles poten- 
tial, one gets 



_Q_ 



2(p+l) 



> 



P {2p + 1) / M 



2tt 
rn 



4+2p 



(68) 



This formula is similar to Eq. (48) of Ref. [34J. It can 
also be re-written as 



(69) 



where T(p) = \p(2p + l)/(2w)] 1 /^ +1 ^ ~ 0(1). Then, 
integrating the Klein-Gordon equation leads to the fol- 
lowing expression 



Qi, 



l-4p(p+l) 

m„, 



l + 2p(p+l) 



M x 


4+2p 

) 


' Qini 








M x 


4+2p 

) 











-2(p+l) 



In 



</>(jV) 



l/2(p+l) 



-2(p+l) 



TV 



l/2(p+l) 



(70) 
(71) 



r 



where in the last equality we have used the fact that the 
evolution of the inflaton field can be approximated by <j> ~ 
0ini[l — N/(2N T )]. The subscript "no inter" just reminds 
that the above equation gives Q in the case where there 
is no interaction between Q and <f>. From this expression, 
we deduce that the quintessence field is frozen if 



A _iQini > x 



(72) 



If this condition is satisfied, then obviously the condi- 
tion (|69|l is also satisfied. The contrary is not necessarily 
true but, as show for instance in Fig. 3 of Ref. Q, this 
only concerns a small range of initial conditions. There- 
fore, we can consider that Q no inter — Qini. 

We are now a position where the values of the 
quintessence field with and without interaction can be 
compared at the end of inflation. Inflation stops when 
the slow-roll parameter e is equal to unity correspond- 
ing to 0end/w pi = l/(2y/n). With the interaction term 
taken into account, the field will be on the attractor Q m in 
and, therefore, its value at the end of inflation is just the 
value of Q m in at the end of inflation, namely 



Qn 



X 



N=N^ 



P 



.167T 2 

= H(p)\. 



l/(2p+2) / 1 



V.20F 



-2/(p+l) 

(73) 
(74) 



As already mentioned, the striking feature of this ex- 
pression is that it does not depend on Qini- For orders of 
magnitude estimate, one can consider that H(p) = 0(1). 
For p = 3, which is our fiducial model, one has Q m i n (N = 



N T ) ~ 1.9 x 10 



-14, 



Therefore the ratio of Q at the 



end of inflation without the interaction term taken into 
account to Q at the end of inflation with the interaction 
term into account is given by 



end 



Qini 



m 
rn 



pi 



v(p+i) 



10 55/(p+l) x >*"'i 



TT \ -1/(P+1) 

^) (75) 



(76) 



This ratio is necessarily greater than one, see Eq. 172fl . 
and for "large" initial conditions can be much bigger than 
one. This means that, generically, Q no intcr ^ Qinton 
i.e. the effect of the interaction is to force the quintessence 
field to remain small during inflation. 

Two loopholes could modify the above conclusion. 
Firstly, we have seen that the numerical integration, on 
which our study is based, is valid only if the initial value 
of Q is not too far from the initial value of the minimum. 
We have used the results obtained under this condition 
and have extrapolated them for any initial conditions. 
If the initial value of the quintessence field is far from 
the time-dependent minimum, we still expect a phase 
of oscillations. However, since we have noticed that the 
amplitude of the oscillations tend to be quite big even if 
the initial displacement from the minimum remains rea- 
sonable (this is not surprising for a potential like Q~ 2p 
which is very "abrupt"), it could actually happen that 
the amplitude of the oscillations, in the case where the 
initial displacement is large, are so big that the assump- 
tion that the quintessence field is a test field becomes 
violated. This is a regime that has not been studied in 
the present article. Some new interesting effects could oc- 
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cur in this case. However, we suspect that the numerical 
study ol this situation could be quite tricky. 

Secondly, it has been shown in Refs. [23| that the quan- 
tum effects could strongly modify the evolution of Q dur- 
ing inflation. These quantum effects have been calculated 
in Refs. (sij by means of the stochastic inflation formal- 
ism for a free field. Therefore, what should now be done 
is to compute these effects in the case where the interac- 
tion term is present. This is clearly a difficult task which 
is beyond the scope of the present paper. In Refs. [34| . 
it has been shown that the quantum effects can push the 
quintessence field to quite large values at the end of in- 
flation. As a consequence, the attractor solution could 
be joined only at late time and even after the present 
time. In this case the quintessential scenario would lose 
an attractive feature, namely its insensitivity to the ini- 
tial conditions. In this respect, the results reached in the 
present article are good news since the effect of the inter- 
action term seems to retain the field to quite small values. 
One could speculate that this could maybe compensate 



the influence of the quantum effects. 

In conclusion we have shown that the non- 
renormalizable interactions between the inflaton and 
the quintessence field have drastic consequences dur- 
ing inflation. The quintessence field follows an attrac- 
tor and remains small compared to the Planck scale 
at the end of inflation. This sets the initial condi- 
tions for the quintessence field. As is well known, the 
quintessence field cannot couple strongly to matter field. 
On the contrary non-renormalizable couplings between 
the quintessence field and matter, in particular Cold 
Dark Matter, may have a crucial impact on the late time 
physics of the quintessence field, i.e. on the coincidence 
problem. This is left for future work. 
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